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ABSTRACT 

Many pulsars are observed to "glitch" , i.e. show sudden jumps in their rotational 
frequency v, some of which can be as large as Au/u « 1CP 6 — 10~ 5 in a subset of pulsars 
known as giant glitchers. Recently Pizzochero (2011) has shown that an analytic model 
based on realistic values for the pinning forces in the crust and for the angular mo- 
mentum transfer in the star can describe the average properties of giant glitches, such 
as the inter-glitch waiting time, the step in frequency and that in frequency derivative. 
In this paper we extend the model (originally developed in Newtonian gravity and for 
a polytropic equation of state) to realistic backgrounds obtained by integrating the 
rclativistic equations of stellar structure and using physically motivated equations of 
state to describe matter in the neutron star. We find that this more detailed treatment 
still reproduces the main features of giant glitches in the Vela pulsar and allows us to 
set constraints on the equation of state. In particular we find that stiffer equations of 
state are favoured and that it is unlikely that the Vela pulsar has a high mass (larger 
than M « 1.5M ). 
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1 INTRODUCTION 

Many pulsars are observed to "glitch" , i.e. they show sudden 
increases in their spin frequency that are instantaneous to 
the accuracy of the data. To date several hundreds of glitches 
have been detected ( jEspinoza et al.|2011| ), with relative in- 
creases in the spin frequency v that range from as low as 
Av/v « 10 -11 to Av/v ~ 10~ 5 . In particular a class of pul- 
sars, of which the Vela pulsar is the prototype, exhibit what 
are known as "giant" glitches ( Espinoza et al.||2011 |, large 
steps in the spin frequency (Av/v m 10 _o ) which are accom- 
panied by an increase in the spindown rate i> and exhibit a 
rough periodicity in their recurrence rate (for example giant 
glitches in the Vela occur roughly every three years). 

Shortly after the first glitches were observed it was sug- 
gested that they could be due to a superfluid component in 
the stellar interior, weakly coupled to the normal component 
and to the electromagnetic emission, that could store angu- 
lar momentum and then release it catastrophically, giving 
rise to a glitch (Baym, Pathick & Pines 1969 Anderson & 
Itoh||1975| |Alpar||1977| |Alpar et al.||1984a| ). Large scale su- 
perfluid components are, in fact, expected in Neutron Star 
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(NS) interiors on theoretical grounds given that the tem- 
perature of the star will drop below the superfluid critical 
temperature (typically « 10 9 K) soon after birth. Further- 
more recent observations of the cooling of the young NS in 
the supernova remnant Cassiopeia A have provided the first 



direct indication of superfluidity in NS interiors ( Page et al 
2011[|Shternin et al.|20TL . 



A superfluid rotates by forming an array of quantised 
vortices which carry the circulation of the fluid. In the NS 
crust the vortices can be strongly attracted, "pinned" , to the 
nuclear lattice ( Anderson fc Itoh|1975 Alpar| 19*77 Pines et 



al.|1980||Anderson et al.|1982[ ) and cannot move outward. If 
the superfluid cannot remove vortices it cannot spin down 
and it therefore acts as an angular momentum reservoir. As 
the crust spins down due to electromagnetic emission a lag 
will develop between the superfluid and the normal compo- 
nent, leading to a hydrodynamical lift force (Magnus force) 
acting on the vortices. Eventually when the lag reaches a 
critical value the pinning force will no longer be able to con- 
trast the hydrodynamical lift and the vortices will unpin, 
transferring their angular momentum to the crust and giv- 
ing rise to a glitch. 

Although there is some evidence that smaller glitches 
in young active pulsars such as the Crab may be related to 
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crust quakes (Crawford & Demianski 2003 
al.|2006| ), there is a growing consensus that the basic picture 



outlined above can be used to describe the main features 
of pulsar glitches. There is still considerable debate on the 
"trigger" for vortex unpinning and several mechanisms have 



been proposed (including crust quakes ( Ruderman 1976 Ru- 
derman, Zhu & Chen||1998|), heat release (Larson & Link 



2002[) and hydrodynamical instabilities ( Glampedakis & An 



dersson||2009| >), but recent work by Warszawski & Melatos 



(2008) (see also Melatos fc Warszawski (2009); Warszawski 



& Melatos (20111) has been very successful in using cellu- 
lar automaton simulations that can track the movement of 
a large number of vortices, to reproduce the distribution of 
glitch sizes and waiting times, and [Haskell, Pizzochero fc| 
Sidery (20121 have produced the first hydrodynamical sim- 
ulation that can follow all stages of a giant glitch, from the 
rise to the relaxation. 

One of the main difficulties in performing such calcu- 
lations, up to now, has been the relative scarcity of real- 
istic estimates of the pinning force between vortices and 
nuclei in the crust. However the recent calculation of IGrilll 



fc Pizzochero (|2012a I (see also Grill (20111; Grill fc Piz 



zochero ( 2012b|>) has filled this gap and produced physically 



consistent pinning profiles that can be used to study pul- 



sar glitches. In fact Pizzochero (20111 has shown that these 



forces can be incorporated into a simple analytical model, 
the so-called "snowplow" model, that predicts the size, step 
in frequency derivative and waiting time of Vela glitches. 
The same forces were used by |Haskell, P izzocher o fc Sidery| 
(2012) n a hydrodynamical model to reproduce also the 
post-glitch relaxation of Vela glitches and show that the 
model is consistent with the size and waiting times of other 
pulsars that exhibit giant glitches. 

We shall discuss the details of the snowplow model in 
the following sections, the main assumption, however, is that 
vortices close to the rotational axis of the star will only be 
weakly pinned at their extremities. The Magnus force can 
thus easily unpin them and move them towards the equator, 
where they will repin as they are now immersed in the strong 
pinning region of the inner crust. This creates a vortex sheet 
that moves close to the maximum of the pinning potential 
as the star spins down. Once the maximum critical lag is 
reached the pinning force can no longer balance the Magnus 
force, and all the vorticity that has been accumulated is free 
to go, giving rise to a giant glitch. This also gives rise to a 
natural periodicity for giant glitches, although crust quakes 
or vortex avalanches can unpin part of the vorticity and give 
rise to smaller glitches before the maximum. 

In this picture we implicitly assume that vortices that 
cross the core of the star are immersed in a low drag environ- 
ment, i.e. that there is no pinning in the core and vortices are 
free to move out. This would not be the case if the protons 
in the interior are in a type II superconducting state and 
there is a strong interaction between magnetic flux tubes 
and rotational vortices ( jRuderman, Zhu &: Chen|1998| |Link| 
20031. Note, however, that a large portion of the star may 
be in a type I superconducting state ( |Jones|[2006[ ) in which 
the magnetic field is not organised in flux tubes, but rather 
in macroscopic regions of normal matter, and the interac- 
tions may be much weaker (Sedrakian 2005) (although see 
Jones| ( |2006[ ) for a discussion of strong interactions in type 
I superconductors). Furthermore recent calculations suggest 



that even if the superconductivity is of type II, the inter- 
action between vortices and flux tubes will be weak in the 
presence of strong entrainment or superfluid E~ hyperons 
| |Babaev||2009[ ). In this paper we thus take the view that 
pinning in the core will be weak, although strong pinning 
of vortices to flux tubes is an intriguing possibility and will 
be the focus of a future publication (Haskell, Pizzochero & 
Seveso, in preparation). 



Pizzochero (20111 and Haskell, Pizzochero & Sidery 



(20121 developed the model for an n = 1 polytrope in New- 



tonian gravity, in order to obtain simple analytical estimates 
for the glitch size. In this paper we investigate the effect on 
the snowplow model of |Pizzochero (20111 of using realis- 



tic Equations of State (EoSs) and relativistic equations for 
the stellar structure. We will see that in general the model 
is consistent with observations and that less massive NSs 
(in the region of 1.2Mq — 1.5Mq) are favoured to describe 
the Vela pulsar. The inclusion of realistic backgrounds and 
equations of state in the hydrodynamical models will be the 
focus of future work. 



2 STELLAR STRUCTURE 



As outlined above the model proposed by Pizzochero (2011 1 



relies on calculating the amount of vorticity that can be 
stored in the inner crust before the Magnus force overcomes 
the maximum pinning force and then calculating the amount 
of angular momentum exchanged between the superfluid 
component and the crust. In order to obtain realistic esti- 
mates of the moment of inertia of the different components 
and of the strength of the density-dependent pinning force 
per unit length, it is thus important to use a realistic equa- 
tion of state (EoS) for dense matter that describes the rela- 
tion P(p) between density and pressure. The stellar struc- 
ture and density profile can then be obtained by integration 
of the Tolman-Oppenheimer-Volkov (TOV) equations: 
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where m(r) is the mass contained in a sphere of radius r, 
p(r) is the density profile and P(r) is the pressure. These 
differential equations model the hydrostatic equilibrium in- 
side the star with relativistic approach and, of course, re- 
quire the P(p(r)) function. The last two equations can be 
combined in one that gives an expression for the mass and 
pressure derivatives and the system can be solved with valid 
initial conditions. We obtain the functions that describe the 
star with the fourth-order Runge-Kutta method, starting 
at r = with m(0) = and p(0) = p c , for a valid choice of 
the central density p c . The integration stops when we reach 
the condition p(R) = 10 _8 p c and we take R as the radius of 
the star. Of course the mass of the star is M = m(R). As 
a result of this integration we have m(r), p(r) and P(r) of 
the star as showed in Fig. [I] 
We use two different EoSs: 



(i) SLy (Douchin & Haensel 20011 is a moderate EoS 
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Figure 1. Mass, pressure and density profiles of 1.4 Mq neu- 
tron star for the EoSs considered: with fixed mass, a stiffer EoS 
produce a star with larger radius and lower central density. 



based on a non-relativistic parametrisation; this equation 
describes the whole star with a single analytical expression 
and so it is more convenient to integrate; 



(ii) GM1 by |Glendenning fc Mos zokowski ( |1991[ ) is a stiff 
P(p) relation that is very similar to SLy in the crust of star, 
but not in the core due to different microscopic approach 
used to describe hadrons at densities higher than po. 

Fig. [2] shows the different mass-central density relations for 
the two equations of state. 

This figure clearly identifies also the presence of a limit 
for the mass for a neutron star. The existence of a maximum 
mass M max is an effect of the relativistic nature of the TOV 
equations, where pressure contributes to the gravitational 
field. Of course M max depends on the EoS: a stiffer equation 
of state gives a higher M max . The evidence of the existence of 
a 2Mq neutron star ( Demorest et al.|2010 | allows, in fact, to 
reject soft equations of state that predict a maximum mass 
below this value. Table[T]shows the different values of M max 
for the equations of state considered. 

In this work we consider stars with masses from IMq 
to M max . With the density profile p(r) it is possible to iden- 
tify, for each star, the structural regions that are relevant 
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Figure 2. This plot shows the mass-central density relation for 
the EoSs considered. As expected, we find a maximum mass value 
above 2Mq for each equation of state (see table III. 



Table 1. This table shows, for each EoS, the maximum allowed 
mass with the corresponding central density p c (in units of nuclear 
saturation density po)> radius of the star R, radius of the core R c 
and radius of the inner crust Ri c . 



EoS 


Pc 


M 


R 


Rc 


Ric 




(po) 


(Mq) 


(km) 


(km) 


(km) 


SLy 


10.2 


2.05 


9.98 


9.68 


9.86 


GM1 


7.1 


2.36 


11.98 


11.57 


11.82 



for the model and important for the pinning. In particu- 
lar we calculate the radius of the core R c as the distance 
from the center of the star where p C orc = p(Rc) ~ 0.5po 
(po = 2.8 x 10 14 gcm -3 is the nuclear saturation density); 
the inner crust-outer crust interface R lc corresponds, on the 
other hand, to the density value pd = 0.0015po that is the 
neutron drip point: this means than in the outer crust there 
are no free neutrons. It is easy to calculate the moment of 
inertia of a shell delimited by radii r\ and T2 : 



I(ri,r 2 ) = 



8tt 



r 4 p(r) dr. 



(4) 



We can then calculate also the moment of inertia of every 
region, considering that /core = 1(0, R c ), lie = I(R c ,Ric) 
and 7 oc = I (Ric, R). 

Table[2]shows all the relevant parameters for the consid- 
ered stars, obtained from the integration of the TOV equa- 
tions with SLy and GM1. 



3 PINNING AND VORTICITY 

One of the most important ingredients of the model is clearly 
/pin(p), the pinning force per unit length that acts on the 
vortex line as a result of its interaction with the lattice (in 
the inner crust). Although the pinning force per pinning site 
can readily be evaluated from the knowledge of the pinning 



energy ( 


Alpar 


1977 


zochero 


2003 


2004 



Epstein & Baym 1988 Donati & Piz- 



2006) 



the force per unit length of a 
vortex, which is the quantity that must be equated to the 
Magnus force in order to understand whether a vortex is 
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Table 2. We give all the structural parameters (as defined in section 2) of the stars used to test the snowplow model, for both 
EoSs tested. See also figure [3] for a graphical representation of these quantities. 



EoS 


M 


R 


Rc 


Ric 


^tot 


I core 


lie 


^OC 






(km) 


(km) 


(km) 


gem ) 


(^lU gem ) 


( i n43 <~c ^w, 2 \ 

(1U gem ) 


ylU gem ) 


SLy 


1.0 


11.86 


10.35 


11.23 


0.739 


0.697 


4.181 


6.638 




1.1 


11.83 


10.49 


11.28 


0.827 


0.788 


3.923 


5.945 




1.2 


11.80 


10.60 


11.31 


0.914 


0.878 


3.652 


5.317 




1.3 


11.76 


10.69 


11.32 


0.999 


0.965 


3.370 


4.738 




1.4 


11.71 


10.75 


11.32 


1.079 


1.048 


3.078 


4.198 




1.5 


11.64 


10.79 


11.29 


1.154 


1.126 


2.777 


3.685 




1.6 


11.55 


10.79 


11.24 


1.222 


1.197 


2.469 


3.194 




1.7 


11.42 


10.76 


11.16 


1.279 


1.258 


2.150 


2.718 




1.8 


11.26 


10.68 


11.03 


1.322 


1.303 


1.818 


2.248 




1.9 


11.03 


10.54 


10.83 


1.339 


1.324 


1.463 


1.769 




2.0 


10.62 


10.23 


10.47 


1.299 


1.289 


1.042 


1.233 


GM1 


1.0 


13.94 


11.79 


13.02 


1.021 


0.896 


12.505 


19.061 




1.1 


13.94 


12.01 


13.12 


1.146 


1.025 


12.068 


17.532 




1.2 


13.94 


12.19 


13.20 


1.271 


1.156 


11.555 


16.108 




1.3 


13.93 


12.35 


13.27 


1.395 


1.285 


10.991 


14.780 




1.4 


13.91 


12.47 


13.32 


1.516 


1.412 


10.382 


13.530 




1.5 


13.89 


12.58 


13.34 


1.634 


1.536 


9.738 


12.340 




1.6 


13.85 


12.66 


13.35 


1.747 


1.657 


9.062 


11.198 




1.7 


13.79 


12.71 


13.35 


1.854 


1.771 


8.362 


10.099 




1.8 


13.72 


12.74 


13.32 


1.954 


1.878 


7.635 


9.031 




1.9 


13.62 


12.74 


13.26 


2.043 


1.974 


6.885 


7.987 




2.0 


13.49 


12.70 


13.17 


2.118 


2.057 


6.107 


6.956 




2.1 


13.33 


12.63 


13.05 


2.173 


2.120 


5.292 


5.922 




2.2 


13.10 


12.48 


12.85 


2.194 


2.150 


4.411 


4.851 




2.3 


12.71 


12.20 


12.51 


2.146 


2.113 


3.371 


3.631 



pinned or free, is much more complex to evaluate, as it de- 
pends on the rigidity of a vortex and on its orientation with 
respect to the crustal lattice. 



Grill fc Pizzochero (2012a I (see also Grill (20111; Grill 



& Pizzochero ( 2012b l) have performed numerical simula- 
tions to evaluate this quantity, taking into account differ- 
ent orientation of the bec lattice. They found that the or- 
der of magnitude of the maximum pinning force /pm is ap- 
proximately 10 1 dyncm -1 and that there is no significant 
difference for the pinning force per unit length in consid- 
ering vortex-nucleus interaction attractive or repulsive in 
different regions. Another interesting result found by |Grill| 
& Pizzochero ( 2012a I regards the position of the maximum 



/pm- In this paper the authors use a density-dependent pair- 
ing gap A(p) obtained with a realistic microscopic nucleon- 
nucleon interaction. It is known that the polarization effects 
of the neutron medium reduce the paring gap, but there is 
yet no agreement on how strong this suppression will be, 
although it seems reasonable to divide the A(p) by a fac- 
tor P between 2 and 3. Grill & Pizzochero (2012a I consider 
the case ft = 1 and /3 = 3 and find that for the two corre- 
sponding profiles /pin(p) the maximum is shifted at differ- 
ent densities, even if the parameter /3 is, of course, only a 
scaling factor on the same pairing gap profile. The precise 
height of the maximum thus depends on the vortex ten- 
sion used in the model (although the order of magnitude 
remains 10 dyncm -1 ) and does not affect the location of 
the maximum (once /3 is fixed). In this work we therefore 
constrain the exact value of the maximum amplitude of the 
pinning force by fitting the average waiting time between 
giant glitches in the Vela pulsar, as will be explained in the 




PtPol 



Figure 4. The profile of the pinning force / p i n (p) for the two 
cases /3 = 1 and /3 = 3, with a choice for the maximum of /pm 
1 dyn cm -1 . 



next sections. In figure [4] we show the two pinning profiles 
/pin(p) used in this work for /3 = 1 and /3 = 3 (plotted here 
with the choice of /pm = 10 15 dyncm -1 ). The case f3i has a 
maximum at p w 0.325po, while in the fiz case the maximum 
is at p ~ 0.14po. In both configurations we take the pinning 
force to vanish at p C ore and pd, due to the fact that the lat- 
tice exists only in the crust and that in the outer crust there 
are no free neutrons to produce vortices. 

A single vortex line will be described parallel to the 
rotational axis and distant from this axis by a distance x, 
that represent the cylindrical radius. We consider also the 
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Figure 3. The first figure shows the dependence of the radius of the neutron star on the total mass, for the SLy and GM1 EoSs. The 
other plots represent the thicknesses of the stellar regions (core, inner crust and outer crust) as function of mass. As one can see a more 
massive star has thinner crusts, while a stiffer equation of state produces a larger star. 



vortex line to be continuous throughout the core: there is, in 
fact, no theoretical evidence for the existence of an interface 
of normal matter between the core and the inner crust, that 
can justify the hypothesis of a core with vorticity separated 



from the crust (Zhou et al.||2004l). Naturally the vortices 



may not be straight and parallel to the rotational axis, as 
turbulence may develop in the stellar interior, especially in 
the presence of strong pinning ( |Lmk||2011b"|a| >. We will not 
consider this possibility here, but will discuss some of its 
likely consequences in the following. 

With the above hypothesis, we can identify two (cylin- 
drical) pinning regions based on the strength of the pinning 
interaction. The strong pinning region is defined by x > R c 
and corresponds to the part of the star in which the vortices 
lie entirely in the inner crust region, and are therefore sub- 
ject to pinning for their whole length. On the other hand, 
in the weak pinning region (x < R c ), a vortex line is pinned 
only at its extremities that are immersed in the crust, while 
there is no pinning interaction in the core (see figure IS|. 



4 THE MODEL 

Thanks to the axial symmetry of the problem, we can de- 
scribe the macroscopic quantities of the rotating superfluid 
in terms of the variable n(x) that represents the number of 
vortices per unit area, at a distance x from the rotational 
axis of the star. The angular velocity Q, s (x) of the super- 



fluid component of the star is in fact proportional to the 
number N(x) of vortices enclosed in a cylindrical region of 
cylindrical radius x, and can be expressed as: 



Z7T X 2 



2lTX 



! n(x') da' 



(5) 



where the integration is performed on the area enclosed by 
the radius x. This result follows from the quantization of the 
circulation per vortex line that is encoded in the constant 
ft = irTi/mN- 

Once a star has been fixed by the choice of an EoS and 
the integration of the TOV equations, the model requires, 
as a first step, the evaluation of the pinning force for the 
whole length of a generic vortex line. This can be obtained 
starting from the function / P m(p) discussed previously. Let 
us imagine a vortex line parallel to the rotational axis of the 
star and distant x: the total pinning acting on it is given by 
the integration of / P i n (p) over its length: 



F pin {x) = 2 



/pin 



•(V 



X 2 + Z' 



dz 



(6) 



where £(x) = I^J R? c — x 2 is the length of the vortex line, 
obtained considering that the vortex line ends at the inner- 
outer crust surface. To understand better the role of the 
pinning force, we choose a neutron star of IAMq with SLy 
equation of state and we plot the function F p i n (x) for x 
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Figure 5. A schematic representation of the geometry of our 
problem (out of scale). The whole shaded area represents the 
inner crust of the NS, where vortices are pinned to the lattice. 
The darker part indicates the strong pinning region, where the 
vortices are subjected to pinning for their whole length. The star 
is threaded by straight continuous vortices. 
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Figure 6. The total pinning force -F p in(:r) integrated on the whole 
length of a vortex line distant x from the rotational axis of the 
star. This plot is obtained taking a star of IAMq and SLy EoS. 
The pinning profile used is the /3 = 1 case plotted in figure [4] 



from to R lc (fig. [6] corresponding to = 1 and /pm = 
lO^dyncm" 1 ). 

The pinning interaction is not the only force that acts on 



a vortex line. As shown in Ruderman & Sutherland ( 19741, 



pinning prevents the vortex line from moving with the local 
superfluid velocity because the vortex line is compelled to 
have the velocity of the normal matter component (the nor- 



mal component rotates as a rigid body with angular velocity 
Q c ). This fact give rise to a Magnus force: 



Kp s e z x (v„ - v s ) 



(7) 



where v„ is the velocity of the vortex line and v 3 is the 
superfluid velocity; here f m must be intended as force per 
unit length. 

In this expression p„ is the density of the superfluid 
fraction of the star. In fact the whole star can be divided 
in two components: the normal one (which includes also the 
protons in the core as they are coupled with the crust by 
the magnetic field) and the superfluid one, on which the 
Magnus force will act. It thus follows that p s = (1 ~ x p )p 
where x p is the proton fraction at a given density. Of course 
this quantity is a microphysical property of matter and for 
this reason is strictly dependent on the EoS used. As this 
information is not provided with the EoSs used, we use the 



results of Zuo et al. ( 2004) who give the proton fraction x v {p) 
as a function of the total density in the case of two-body 
interactions and also in the case of three-body forces. We use 



both the x p (p) relations of Zuo et al. (20041 but we consider 
also a third case where the proton fraction is a constant that 
does not depend on the total density. We also introduce the 
parameter Q that represent the superfluid fraction of the 
star. It is defined for the general case as 



Q = Is/Iu 



!o r 4 p{r) dr 



(8) 



where we have used eq. |4j I to t is the total moment of inertia 
and I s is the moment of inertia of the superfluid component. 
In the case of a constant proton fraction it then follows that 
Q = 1 — x p . The average value is Q ~ 0.95 and therefore we 
shall test our model also with this prescription. 

The Magnus force in [7] has only one component in the 
radial direction (v„ and v s are, in fact, directed along eg, 
so the cross product is directed along e x ) and therefore can 
be rewritten as: 



t m {x, z) = f m (x, z) e x = -np s (x, z)xAfl(x) e x 



(9) 



where the difference of the two velocities is written as 
Av(x) = xAfl(x) = x [£l c — £l s (x)] and depends only on 
the coordinate x, as described by equation [5] This quantity 
is negative between two glitches because the normal compo- 
nent spins slower that the superfluid one; indeed the Magnus 
force is a hydrodynamical lift that pushes the vortex outward 
from the rotational axis. The key point here is the fact that 
the normal component spins down as a consequence of the 
loss of energy by electromagnetic radiation of the star: the 
result is an increase of fm(x, z) in the time interval between 
glitches. 

The same integration performed with / p i n over the 
length of the vortex can be done with the Magnus force. 
We can consider the total Magnus force acting on a vortex 
line distant x from the rotational axis: 



F m {x) = 2 



l(x)/2 



fm(x, z) d,Z 



= 2kxA£1{x) 



Ufx 2 



dz. 



(10) 



The basic idea here is to compare the pinning force and 
the Magnus force to find the critical lag Af2 C r(s) that repre- 
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1.2-10 24 




forces and for the short relaxation timescale r c ~ 1 — 10s 



Figure 7. Plot of the expression F^(x) = F m (x)/ AQ(x) as a 
function of the cylindrical radius x. This plot is obtained taking 
a star of 1.4M© and SLy EoS. 



sents the depinning condition: when the actual lag between 
the two components of the stars reaches the value Afi cr at 
some point with cylindrical radius x. the vortices here are 
unbound from the lattice, as the Magnus force now exceeds 
the pinning interaction that held the vortices in place: 

F pin (x) = F m (x) = ACl CI (x)F^(x) (11) 

Here F^(x) = F m (x) / AQ. cr (x) and it is plotted in fig. [7] 
using the same reference star as in fig. [6] The important 
quantity is therefore the critical lag that can be easily eval- 
uated as: 



AQ CI (x) = 



KX Jo ps W x 



dz 



(12) 



p s I Va^ + z 2 ) dz 

In figure [5] we plot the critical lag for sample stars from 
table O 

It is important to point out that the critical lag shows 
a peak Af2 crmax = Afi cr (a^ max ) in a region that corresponds 
to the inner crust, that is the region where the pinning is 
stronger. In this region our estimate of Af2 cr (a;) is reasonable 
since pinning is continuous along the whole single vortex. 
This is not the case for the critical lag in the core, because 
here pinning acts on vortices only at their extremities: as 



explained by Pizzochero ( 2011 1 this fact is responsible of the 



weak pinning in this region, even though we can assume that 
the system maintains axial symmetry due to the collective 
rigidity of vortex bundles ( Ruderman fc Sutherland) 1974 1. 

As the star slows down, the depinning condition 
Af2(x) ^ AQ cr (x) is first reached in the core: as shown byfl 



Link ( 2009 1 , in this region repinning is dynamically possible 



if the lag falls below a critical value (smaller than the one for 
depinning). This suggest the following interpretation: in the 
core, as the star slows down, the vortices are continuously 
depinned and repinned, establishing a dynamical creep that 
removes the excess vorticity on short timescales. Further- 
more the Magnus force in the interior is likely to overcome 
the tension of vortices and depin them long before the un- 
pinning conditio n in the crust is met (|Adams, C ieplak & 



Glaberson 1984 Haskell, Pizzochero & Sidery 20121. The 



conclusion is thus that vortices in the core can essentially be 
considered free. In this region the scattering of electrons off 
magnetised vortex cores is mainly responsible for the drag 



( Alpar et al.|1984b||Andersson, Sidery fc Comer|2006[ ): this 
means that we can consider the normal and the superffuid 
components in the core as coupled with a lag of order |f2|r c . 

In the time between glitches, the depinning region be- 
comes larger, involving also the crust: in the inner crust the 
excess vorticity is repinned and creates a thin vortex sheet 
that moves toward the peak: this sheet is pushed outward by 
the increasing Magnus force and it stores angular momen- 
tum. When the peak is reached, there is no more pinning 
interaction that can block the excess vorticity: this vortic- 
ity is suddenly released and reaches the outer crust. At this 
moment the angular momentum stored by vortices is trans- 
ferred to the normal component of the star, and this causes 
the glitch. 

It is straightforward now to evaluate the time interval 
between two glitches: this is given by the time needed to 
create a lag Af2 crmax 



At s 



AQ CI 



\Q\ 



(13) 



where Q is the deceleration of the normal component re- 
ferred to the pre-glitch steady-state condition. 

The above arguments indicates that, immediately be- 
fore a glitch, a lag of Af2 crmax will have been created for 
x = i max . This means that we can use equation [5] to express 
the number of vortices stored at the peak in the sheet just 
before a glitch: 



2tt 



(14) 



Due to the particular shape of the critical lag in figure [8j we 
can assume that in this moment the excess vorticity in the 
region x > x max has been entirely removed by the Magnus 
force, and therefore the N v vortices are the only ones respon- 
sible for the transfer of angular momentum to the normal 
component of the star. To evaluate the angular momentum 
transfer we start from the definition dL = fl 3 (x) dl s . As we 
are interested in the angular momentum stored by N v vor- 
tices at the peak of the pinning potential, we use the relation 
in |14| and perform the integration on the cylindrical region 
Xmax < x < R c to obtain the requested quantity (the inte- 
gration on the coordinate x stops at R c due to the fact that 
in the outer crust there is no superfluid component): 



AL g i = 2kN v 



Rc 



x dx 



x 2 + z 2 ) dz. 



(15) 



Following the arguments above, at the moment of a 
glitch only a fraction of the core superfluidity is coupled 
to the normal component of the star: in fact, the rise time 
of a glitch (r g i) is very short and only the instantaneously 
depinned fraction of vorticity in the core can respond to the 
variation of the angular velocity of the crust. We introduce 
the parameter Y g i to encode this fractional quantity. In fact 
the best observational upper limit is of r g i < 40 s for the 
Vela 2000 glitch ( podson, McCulloch fc Lewis||2002[ ), while 
an interesting lower limit of r g i > 10~ ms can be set by 
the non-detection of gravitational waves from the Vela 2006 
glitch (Warszawski & Melatos 20121. Theoretical estimates 



suggest that r g i 
20121. 
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Figure 8. Plot of the critical lag AQ CT (x) for different stellar models, with varying mass and equations of state. The pinning profiles 
used are plotted in figure [4] and we consider both the case f) = 1 and the case f) = 3. Note that in both cases ((3 = 1 and /3 = 3) the 
maximum amplitude of the pinning force is the same, so the difference in the maximum lag for the two cases is now entirely due to the 
different position in density of the maximum in the pinning force profile. 



In the pre-glitch steady-state condition, due to the long 
timescales involved, we can assume Yoo = 1; but during a 
glitch this quantity cannot be calculated with the snowplow 
model as it depends on the detailed short-time dynamics of 
the vortices, and must thus be determined with hydrody- 
namical simulations such as those in [Haskell, Pizzoc hero fc| 
Sidery ( 2012 1. As this is beyond the scope of the current pa- 



per, the quantity Y g \ is taken as a parameter of this model, 
and must be inferred from the observational data as shown 
in the next section. 

The value Y R \ is needed for the evaluation of Af2 g i, the 
jump in angular velocity of the normal component of the star 



due to a glitch. This corresponds to the ratio between the 
angular momentum transfer AL g i and the effective moment 
of inertia I e g of the coupled fraction of matter during the 
glitch. One thus has that I e g = (1 — Q)I to t + QY g \I to t and 
the requested quantity is therefore: 



AfLi = 



AL Rl 



/tot [1 - Q (1 - Y gl )] ■ 



(16) 



A further parameter of the glitch that can be calculated 
is the relative acceleration of the crust. As illustrated in 



Pizzochero (2011 1 the desired relation follows from variation 
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at the glitch of the Euler equation for the normal component 
and angular momentum conservation: 



AQg 



3(1 - Y sl ) 
1 - Q(l - Vg,) 



(17) 



5 RESULTS AND OBSERVATIONS 

In this section we test the model proposed here against 
observations. As the model has been developed for giant 
glitches we shall compare our results to observations of gi- 
ant glitches in the Vela pulsar. The Vela (PSR B0833-45 
or PSR J0835-4510) has a spin frequency v w 11.19 Hz 
and spin-down rate v ~ —1.55 x 1CP 11 Hz s _1 ; from re- 
lation [13] this value correspond to a maximum critical lag 
of Ar2 crmax = 8.6 x 10 -3 rads -1 , where we have consid- 
ered that the average time between glitches for this pul- 
sar is 2.8 years. The glitch is usually described in terms 
of permanent steps in the frequency and frequency deriva- 
tive and a series of transient terms that decay exponen- 
tially. It is well known that at least three transient terms 
are required, with decay timescales that range from months 
to hours (Flanagan 1996). Recent observations of the 2000 
and 2004 glitch have shown that an additional term is re- 
quired on short timescales, with a decay time of approx- 
imately a minute. Given that the detection in 2004 was 
only barely above the noise we shall refer to the January 
2000 glitch. In this case the j ump in angular velocity was 
of Af2 g i = 2.2 x 10~ 4 rad s" 1 ( |Dodson, McCulloch fc Lewis 
2002||Dodson et al.||2007[ ). This is a fairly typical value for 
giant glitches in the Vela, and we take it as our reference 
value. The relative step in frequency derivative correspond- 
ing to the transient term with the shortest decay timescale 
(1 minute) for this glitch is AQgi/Q^ ~ 18 ± 6 (la error), 
and we assume that this is a good approximation to the in- 
stantaneous step in frequency derivative at the time of the 
glitch. 

As explained in the previous section, the model has two 
free parameters that are the maximum value of the pin- 
ning force value /pm and Y g \: this means that, once a star 
has been fixed (by choosing the EoS, the mass M, and the 
superfluid fraction relation) we can use two observational 
quantities to constrain the parameters of the model and 
compare further observables to the quantities predicted by 
calculations. In particular, for each fixed star, we rescale the 
maximum of the pinning force in order to produce the maxi- 
mum critical lag Af2 crmax required to reproduce the average 
waiting time between glitches in the Vela. This allows us 
to calculate directly and univocally the angular momentum 
AL g i from equation [15] As we want to reprodu ce a glitch 
of amplitude Af2 g i = 2.2 x 10 -4 rads -1 , relation 
rewritten in the following form 



16 



can be 



AL„ 



/tot Ai^p 



+ 1 



(18) 



and therefore can be used to fix the coupled fraction of mat- 
ter during the glitch. Tables[3][4]and[5]give the fitting param- 
eters for all the configurations tested. We can see that the 
value of the maximum pinning force /pm does not change 
significantly with the total mass of the star. In these ta- 
bles negative values of Y s \ are not given as they would not 



be physically acceptable: a negative value would mean that 
there is not enough angular momentum to produce the re- 
quired jump in angular velocity, even if we consider the core 
vorticity completely decoupled from the normal component 
of the star at the time of the glitch. 

The remaining tables, numbered [6] [7] and [8] show the 
physical quantities that the "snowplow" model permits to 
evaluate. These are of course the angular momentum AL g i 
transferred to the crust during the glitch and the relative 
step in frequency derivative. We can see that the order of 
magnitude for AL g i is 10 40 ergs, that is compatible with 
the upper limits on the glitch energy obtained from obser- 



vations of the power wind nebula surrounding Vela ( Helfand 
et al.|200T l and with the results found in Pizzochero (2011 1, 



where the same model is applied analytically with a poly- 
tropic EoS in Newtonian gravity. From these tables one can 
see that, for a particular choice of EoS and proton fraction, 
the angular momentum AL g i stored by vortices decreases 
with the total mass of the star. This behaviour can be easily 



explained, as shown by Pizzochero (20111, in terms of the 
quantity £ max /i?i c shown in the tables: AL g i is obviously 
related to the number N v of vortices stored at the peak (in 
tables [6] [7] and [8] see also eq. 1 1 5b - that however doesn't 
change significantly with the mass - but it depends strongly 
on the ratio £max/-Ric which increases at higher masses. In 



Pizzochero (20111 (fig. 4) it is clearly shown that the angu- 



lar momentum stored by the vortices at the peak decreases 
rapidly moving the position of the peak towards the outer 
crust. The quantity AL g i also depends on the equation of 
state used (a stiffer EoS produces higher values of AL g i) 
and on the pinning profile: the j3 — 3 condition, when other 
variables are fixed, gives lower values for the angular mo- 
mentum, accordingly to the fact that the relative position 
of the peak with respect to the inner crust radius is higher. 

The "snowplow" model permits to calculate also the 
step in spin-down rate immediately after a glitch, and this 
quantity is given in our tables as AQgi/Cloo. It has been 
calculated only for acceptable values of Y s i, and must be 
compared with the reference value of AQgi/Q^ = 18 ± 6, 



taken from the Vela 2000 glitch ( |Dodson, McCulloch fe 
Lewis|20"02 I. These values suggest that the /3 = 3 configura- 



tions are preferred and this can be considered in reasonable 
agreement with the microscopic results found by |Gandolfi| 



et al. (20081: they find that a realistic suppression factor for 
the pairing gap A(p) is /3 ~ 1.5 but, crucially, also that the 
maximum for A(p) is shifted at lower densities. This leads 
to a profile close to what we obtain for /3 = 3 in our model. 



Finally let us remark that the results in tables [6] [7] and 
[8] for the (microscopically favoured) case j3 — 3, seem to 
indicate that a stiffer equation of state (GM1) is preferred 
as is a lower mass (possibly in the region of IAMq) for 
the Vela pulsar. Naturally such a quantitative conclusion 
is difficult to make on the basis of one observation and it 
would be highly desirable to have information on the short- 
timescale post-glitch behaviour not only of other Vela giant 
glitches, but also of other glitching pulsars. Note that short 
term components of the relaxation have not been measured 
for other giant glitchers, however the "snowplow" model can 
be used to predict waiting times, obtaining results which are 
consistent with observations ( jHaskell, Pizzochero fc Sidery| 
2012 1. 
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Table 3. This table gives the fitting parameters /pm (maximum of the pinning force per unit length) and Y g \ (fraction of coupled 
vorticity at the glitch) for all the considered configurations. a; max is the position (cylindrical radius) of the maximum critical lag. 
These values refer to a constant proton fraction x p (p) = 0.05 (the corresponding fraction of moment of inertia due to the superfluid 
component of the star is Q = 0.95). Unphysical (negative) values for Y g \ are not reported (see text for details). Since the angular 
momentum transferred to the crust during a glitch is strongly dependant on the ratio x max /fiic (see section Hp, this quantity is 
also reported in table. 



13 = 1 /9 = 3 



EoS 


M 
(Ms) 


3^ max 

(km) 


^max / ^ic 


/pm/10 15 
(dyncm -1 ) 




^max 

(km) 


3<max/ ^?ic 


/pm/10 15 
(dyncm - 1 ) 




SLv 


1.0 


10.530 


0.938 


1.562 


0.203 


10.724 


0.955 


0.697 


0.027 




1.1 


10.654 


0.945 


1.581 


0.148 


10.829 


0.960 


0.704 


0.009 




1.2 


10.753 


0.951 


1.595 


0.106 


10.910 


0.965 


0.709 






1.3 


10.827 


0.956 


1.606 


0.074 


10.968 


0.968 


0.713 






1. 1 


10.875 


0.961 


1.613 


0.048 


11.001 


0.972 


0.715 






1.5 


10.897 


0.965 


1.616 


0.027 


11.010 


0.975 


0.716 






1 .u 


1 


u.yuy 




n m 1 
u.u r 1 


1 n QQn 
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n Q77 








1.7 


10.847 


0.972 


1.609 




10.937 


0.980 


0.711 






1.8 


10.759 


0.975 


1.596 




10.838 


0.982 


0.705 






1.9 


10.600 


0.979 


1.572 




10.667 


0.985 


0.694 






2.0 


10.279 


0.982 


1.525 




10.332 


0.987 


0.672 




GM1 


1.0 


12.129 


0.932 


1.798 


0.493 


12.447 


0.956 


0.809 


0.078 




1.1 


12.315 


0.939 


1.825 


0.389 


12.604 


0.961 


0.819 


0.051 




1.2 


12.473 


0.945 


1.849 


0.307 


12.737 


0.965 


0.827 


0.031 




1.3 


12.604 


0.950 


1.868 


0.242 


12.844 


0.968 


0.834 


0.015 




1.1 


12.710 


0.955 


1.884 


0.190 


12.929 


0.971 


0.840 


0.003 




1.5 


12.792 


0.959 


1.896 


0.148 


12.992 


0.974 


0.844 






1.6 


12.852 


0.962 


1.905 


0.113 


13.034 


0.976 


0.847 






1.7 


12.890 


0.966 


1.910 


0.083 


13.055 


0.978 


0.848 






1.8 


12.901 


0.969 


1.912 


0.060 


13.052 


0.980 


0.848 






1.9 


12.885 


0.972 


1.910 


0.040 


13.022 


0.982 


0.846 






2.0 


12.836 


0.974 


1.902 


0.023 


12.960 


0.984 


0.842 






2.1 


12.744 


0.977 


1.889 


0.007 


12.854 


0.985 


0.835 






2.2 


12.586 


0.980 


1.865 




12.683 


0.987 


0.824 






2. .3 


12.287 


0.982 


1.821 




12.367 


0.989 


0.803 





6 CONCLUSIONS 



In this paper we have extended the "snowplow" model for 



giant pulsar glitches of Pizzochero (2011 1 to incorporate rel- 



ativists background stellar models and realistic equations of 
state. In particular we test the model for the SLy and GM1 
equations of state. Unfortunately these equations of state do 
not include information on beta equilibrium, so we use the 
proton fractions calculated by |Zuo et ah] ( |2004[ ) . It would 
of course be highly desirable to use proton fractions that 
are consistent with the individual equations of state in fu- 
ture work, in order to set stringent constraints. Furthermore 
we use, for the first time, the realistic profiles for the pin- 
ning forces per unit length cal culated by [Grill fc Pizzochero| 
( |2012a| ) (see also [GrllT] ( poTT] ); | Grill fc Pizzochero| ( |2012b[ )), 
in order to evaluate the amount of angular momentum that 
can be transferred to the crust during a glitch. 

The model contains three free parameters, the mass of 
the star M, the fraction of superfluid that is coupled to 
the crust during a glitch Y g i , (which can only be estimated 
with dynamical simulations such as those of |Haskell, Piz-| 
|zochero fc Sidery] ( |2012| )) , and the maximum amplitude of 
the pinning force, /pm- Note in fact that while the location 
of the maximum is precisely determined by the microphysi- 
cal calc ulations of [Grill fe Pizzochero] ( |2012a[ ) (see also [Grill] 
(2011|);|Grill & Pizzochero| (|2012b|), the actual value of the 



maximum can vary by factors of order unity or more as it 
depends on the poorly constrained value of the vortex ten- 
sion. We thus treat it as a normalization and determine its 
value by requiring that the waiting time between glitches is 
of 2.8 years, as is approximately the case for Vela glitches. 
We then fit the size of the glitch to an average Vela glitch 
to obtain the value of Y g i. In particular we take the value 
of the Vela 2000 glitch, A » = 2.2 x 10" 4 rads" 1 ( |Dodson, 
McCulloch fc Lewis|2002 i. 



Having determined the free parameters in our model, 
except for the mass of the NS which is free, we compare our 
results to the post glitch step in frequency derivative. Unfor- 
tunately the changes in v on short time scales after a glitch 
are observationally challenging to detect and it has been 
possible to fit for transient steps in frequency and frequency 
derivative on timescales of minutes after a glitch only for 



the Vela 2000 and 2004 gl itch ( podson, McCulloch fc Lewis 
|2002[ |Dodson et al.||2007|j. Given that the detection is only 



barely above the noise for the 2004 glitch ( |Dodson et al. 



2007) we fit to the values obtained for the 2000 glitch, which 



we assume to be a good approximation of the instantaneous 
post glitch behaviour. This justifies our choice of also fitting 
to the value of the jump in frequency of the Vela 2000 glitch. 

The comparison of the model to the observational con- 
straints first of all highlights that the general results of the 



analytic model of Pizzochero (20111 remain valid even in 



Realistic EOSs and the "snowplow" glitch model. 



11 



Table 4. The fitting parameters given here (defined in table |3J| refer to the proton fraction x p (p) proposed by |Zuo et al.||2004| 
obtained with two-body forces. 
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0.945 
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10.753 
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10.910 
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_ 
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10.827 
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_ 
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0.746 


_ 
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1.673 


_ 
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0.745 
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_ 


10.937 


0.980 


0.742 


_ 
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10.759 
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_ 


10.837 
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_ 




1.9 


10.600 


0.979 


1.628 


_ 


10.666 


0.985 


0.723 






2.0 


10.279 


0.982 


1.579 




10.332 


0.987 


0.701 




GM1 


1.0 


12.129 


0.932 


1.861 


0.522 


12.446 


0.956 


0.843 


0.099 




1.1 


12.315 


0.939 


1.890 


0.415 


12.605 


0.961 


0.854 


0.070 




1.2 


12.473 


0.945 


1.914 


0.330 


12.736 


0.965 


0.863 


0.048 




1.3 


12.604 


0.950 


1.934 


0.262 


12.844 


0.968 


0.870 


0.029 




1.4 


12.710 


0.955 


1.951 


0.207 


12.929 


0.971 


0.876 


0.014 




1.5 


12.792 


0.959 


1.963 


0.162 


12.992 


0.974 


0.880 


0.002 




1.6 


12.852 


0.962 


1.972 


0.124 


13.034 


0.976 


0.883 






1.7 


12.888 


0.966 


1.978 


0.092 


13.055 


0.978 


0.884 






1.8 


12.900 


0.969 


1.980 


0.064 


13.052 


0.980 


0.884 






1.9 


12.885 


0.972 


1.977 


0.039 


13.022 


0.982 


0.882 






2.0 


12.836 


0.974 


1.970 


0.018 


12.960 


0.984 


0.878 






2.1 


12.744 


0.977 


1.956 




12.854 


0.985 


0.871 






2.2 


12.586 


0.980 


1.931 




12.683 


0.987 


0.859 






2.. 3 


12.286 


0.982 


1.885 




12.367 


0.989 


0.838 





our more physically realistic approach and the results are 
in general consistent for both equations of state for a rea- 
sonable range of neutron star masses. The glitch model pre- 
sented here thus appears robust and compatible with the 
observations of giant glitches in the Vela and is, as shown 
in 



star appears to be justified by microphysical estimates, that 
do not predict an interface of normal matter between the 
crust and core superfluid ( Zhou et al.p 004), it may be the 
case that if the protons in the core are in a type II supercon- 
ducting state this could lead to strong pinning also in the 



Haskell, Pizzochero & Sidery (20121, compatible with stellar interior (Ruderman, Zhu & Chen 1998 Link 20031 



the average waiting time between giant glitches in other 
pulsars. This further reinforces the hypothesis that giant 
glitches are approximately periodic phenomena that occur 
close to the maximum lag that the pinning force can sup- 
port in the crust, while smaller glitches may be triggered by 



Melatos 2011 



random events such as crust quakes ( Ruderman|l976| |Rud 
erman, Zhu fe Chen|1998 \ or vortex avalanch es (jWarszawski 



fc Melatos| 2008 Me latos fc Warszawski|2009||Warszawski fc 



2012|. Furthermore our results favour lower 



masses for the Vela pulsar (smaller than 1.5Mq) and stiffer 
equations of state. Note however that such a quantitative 
conclusion is difficult to draw as not only are we comparing 
to a single observation but dynamical simulations have also 
shown that superfluid mutual friction will contribute sig- 



nificantly to the short term post-glitch spindown ( Haskell, 
|Pizzochero fc Sidery|20T2| as may friction at the crust/core 
interface (van Eysden & Melatos 2010). In will thus be nec- 
essary to further develop hydrodynamical glitch simulations 
to truly constrain the equation of state and the stellar mass. 

Finally let us remark that in this paper we have as- 
sumed straight vortices that cross cross the core of the neu- 
tron star and are only weakly pinned at their extremities. 
Although the assumption of vortices that pass through the 



In this case not only would vortex motion be impeded, but 
it is also likely that turbulence will develop ( |Lhik|2011b|a |. 
Note, however, that a large portion of the star may be in 



a type I superconducting state ( Jones 2006 1 in which the 



magnetic field is not organised in flux tubes, but rather in 
macroscopic regions of normal matter, and the interactions 



may be much weaker ( Sedrakian|2005 1 (although see Jones 



( 2006 1 for a discussion of pinning in type I superconduc- 



tors). Furthermore Babaev (20091 has recently shown that 
in the presence of strong entrainment or superfluid E~ hy- 
perons, the interaction between vortices and flux tubes may 
be much weaker than generally assumed in the presence of 
type II superconductivity. In this paper we thus take the 
view that pinning in the core will be weak, although strong 
pinning of vortices to flux tubes is an intriguing possibil- 
ity and will be the focus of a future publication (Haskell, 
Pizzochero & Seveso, in preparation). 

Turbulence, on the other hand, is well known from labo- 
ratory superfluids and may play an important role in pulsar 



glitches (Peralta et al. 2006 Melatos & Peralta 2007 Per 



alta & Melatos 2009 ) and could couple the superfluid and the 



normal component on inter-glitch timescales in the presence 
of core pinning ( Andersson, Sidery fc Comer|200T l. The in- 
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Table 5. This table is analogous to tables [3] and |4| here the proton fraction used is that calculated by |Zuo et al.| j2004} with 
three-body forces. 



(5 = 1 



f} = 3 



XT' O 


ivi 


% max 


1- 1 f? 

-^max/ -fi-ic 


/PM/1U 


r gl 


% max 


~, / D 

■'■max/ Ji-ic 


7PM/ 10 






( M e) 


(km ) 




(dyn cm~^) 




(km) 




(dyn cm -1 ) 




SLy 


1.0 


10.530 


0.938 


1.619 


0.197 


10.724 


0.955 


0.728 


0.011 




1.1 


10.654 


0.945 


1.638 


0.133 


10.829 


0.960 


0.735 


_ 




1.2 


10.753 


0.951 


1.653 


0.080 


10.910 


0.965 


0.740 


_ 




1.3 


10.827 


0.956 


1.664 


0.036 


10.968 


0.968 


0.744 


_ 




i.i 


10.875 


0.961 


1.672 


_ 


11.001 


0.972 


0.746 


_ 




1.5 


10.897 


0.965 


1.675 




11.010 


0.975 


0.747 


_ 




1.6 


10.889 


0.969 


1.674 


_ 


10.990 


0.977 


0.746 


_ 




1.7 


10.847 


0.972 


1.667 


_ 


10.937 


0.980 


0.742 


_ 




1.8 


10.759 


0.975 


1.654 


_ 


10.837 


0.982 


0.735 


_ 




1.9 


10.600 


0.979 


1.629 


_ 


10.666 


0.985 


0.724 


_ 




2.0 


10.279 


0.982 


1.580 




10.332 


0.987 


0.701 




GM1 


1.0 


12.129 


0.932 


1.862 


0.521 


12.446 


0.956 


0.843 


0.095 




1.1 


12.315 


0.939 


1.891 


0.413 


12.605 


0.961 


0.854 


0.065 




1.2 


12.473 


0.945 


1.915 


0.326 


12.736 


0.965 


0.863 


0.042 




1.3 


12.604 


0.950 


1.935 


0.257 


12.844 


0.968 


0.870 


0.022 




1.1 


12.710 


0.955 


1.952 


0.200 


12.929 


0.971 


0.876 


0.006 




1.5 


12.792 


0.959 


1.964 


0.153 


12.992 


0.974 


0.880 






1.6 


12.852 


0.962 


1.973 


0.112 


13.034 


0.976 


0.883 






1.7 


12.888 


0.966 


1.979 


0.077 


13.055 


0.978 


0.885 






1.8 


12.900 


0.969 


1.981 


0.046 


13.052 


0.980 


0.884 






1.9 


12.885 


0.972 


1.978 


0.017 


13.022 


0.982 


0.882 






2.0 


12.836 


0.974 


1.971 




12.960 


0.984 


0.878 






2.1 


12.744 


0.977 


1.957 




12.854 


0.985 


0.871 






2.2 


12.586 


0.980 


1.932 




12.683 


0.987 


0.859 






2.. 3 


12.286 


0.982 


1.886 




12.367 


0.989 


0.838 





elusion of turbulence in a hydrodynamical glitch simulations 
is, however, a complex matter as not only is the nature of the 



turbulence not known (see e.g. Andersson, Sidery & Comer 
j 2007| |; [Link] ( |2011b|a| )) but also the definition of winning 
force per unit length must be revisited in the presence of 
a turbulent tangle. Such a fundamental issue should clearly 
be the focus of future work. 
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Table 7. The quantities AL g i and AQgi/Cloo here reported follows from calculation base on the proton fraction proposed by Zuo 
|et al.| |2004[ ) with two— body interactions. The corresponding fitting parameters are shown in table [4] Q = I s /Itot is the global 
supcrfluid fraction of moment of inertia. 



fS = l 



EoS 


M 




R 


Q 


N v 


AL gl 




A'v 


AL gl 








(km) 




(10 13 ) 


(10 40 ergs) 


(10 13 ) 


(10 40 ergs) 


SLy 


1.0 


11 


855 


0.948 


3.041 


4.048 


2.962 


3.154 


1.272 


11.609 




1.1 


11 


830 


0.945 


3.113 


3.565 


4.036 


3.216 


1.108 


15.210 




1.2 


11 


797 


0.942 


3.171 


3.119 


5.362 


3.265 


0.961 






1.3 


11 


758 


0.938 


3.215 


2.708 


7.005 


3.299 


0.828 






1.4 


11 


705 


0.934 


3.244 


2.327 


9.065 


3.319 


0.708 






1.5 


11 


635 


0.930 


3.257 


1.981 


11.650 


3.325 


0.596 






1.6 


11 


545 


0.925 


3.252 


1.660 




3.313 


0.499 






1.7 


11 


422 


0.920 


3.227 


1.357 




3.280 


0.405 






1.8 


11 


260 


0.913 


3.175 


1.076 




3.221 


0.321 






1.9 


11 


025 


0.905 


3.082 


0.803 




3.120 


0.239 






2.0 


10 


620 


0.890 


2.898 


0.518 




2.928 


0.153 




GM1 


1.0 


13 


940 


0.966 


4.034 


11.939 


0.856 


4.248 


2.874 


6.713 




1.1 


13 


943 


0.964 


4.159 


10.850 


1.292 


4.357 


2.559 


8.719 




1.2 


13 


940 


0.962 


4.267 


9.805 


1.814 


4.448 


2.303 


10.982 




1.3 


13 


930 


0.961 


4.357 


8.820 


2.433 


4.524 


2.047 


13.792 




1.4 


13 


913 


0.959 


4.430 


7.903 


3.164 


4.584 


1.822 


17.063 




1.5 


13 


885 


0.956 


4.488 


7.034 


4.041 


4.629 


1.612 


20.994 




1.6 


13 


845 


0.954 


4.530 


6.216 


5.101 


4.660 


1.416 






1.7 


13 


788 


0.952 


4.556 


5.463 


6.367 


4.674 


1.237 






1.8 


13 


715 


0.949 


4.564 


4.743 


7.943 


4.672 


1.066 






1.9 


13 


620 


0.946 


4.553 


4.051 


9.948 


4.651 


0.909 






2.0 


13 


495 


0.942 


4.519 


3.415 


12.463 


4.606 


0.760 






2.1 


13 


330 


0.938 


4.455 


2.796 




4.532 


0.620 






2.2 


13 


095 


0.933 


4.345 


2.188 




4.411 


0.481 






2.3 


12 


713 


0.925 


4.140 


1.552 




4.195 


0.342 





Realistic EOSs and the "snowplow" glitch model. 



15 



Table 8. The proton fraction used for this table is the three-body forces model of |Zuo et al.| ( p004} . The corresponding fitting 
parameters are shown in table [5] 
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0.950 


4.430 


7.906 


3.163 


4.584 


1.822 


17.061 




1.5 


13 


.885 


0.946 


4.488 


7.037 


4.039 


4.629 


1.612 






1.6 


13 


.845 


0.942 


4.530 


6.218 


5.099 


4.660 


1.416 






1.7 


13 


.788 


0.937 


4.556 


5.465 


6.365 


4.674 


1.237 






1.8 


13 


.715 


0.931 


4.564 


4.744 


7.939 


4.672 


1.066 






1.9 


13 


.620 


0.925 


4.553 


4.052 


9.944 


4.651 


0.909 






2.0 


13 


.495 


0.917 


4.519 


3.416 




4.606 


0.760 






2.1 


13 


.330 


0.907 


4.455 


2.797 




4.532 


0.620 






2.2 


13 


.095 


0.894 


4.345 


2.188 




4.411 


0.481 






2.3 


12 


.713 


0.873 


4.140 


1.553 




4.195 


0.342 





